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The Faraday problem is an important pattern-forming system that provides some middle 
ground between systems where the initial instability involves just a single mode and in 
which complexity then results from mode interactions or secondary bifurcations, and 
cases where a system is highly turbulent and many spatial and temporal modes are 
excited. It has been a rich source of novel patterns and of theoretical work aimed at 
understanding how and why such patterns occur. Yet it is particularly challenging to tie 
theory to experiment: the experiments are difficult to perform; the parameter regime of 
interest (large box, moderate viscosity) along with the technical difficulties of solving the 
free boundary Navier-Stokes equations make numerical solution of the problem hard; and 
the fact that the instabilities result in an entire circle of unstable wavevectors presents 
considerable theoretical difficulties. 

In principle, weakly nonlinear theory should be able to predict which patterns are stable 
near pattern onset. In this paper we present the first quantitative comparison between 
weakly nonlinear theory of the full Navier-Stokes equations and (previously published) 
experimental results for the Faraday problem with multiple frequency forcing. We confirm 
that three-wave interactions sit at the heart of why complex patterns are stablised but 
also highlight some discrepancies between theory and experiment. These suggest the need 
for further experimental and theoretical work to fully investigate the issues of pattern 
bistability and the role of bicritical/tricritical points in determining bifurcation structure. 
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1. Introduction 


Since Faraday (1831) identified that regular patterns can appear on the surface of a 


shaken container of fluid, the Faraday experiment has been an important system for 
investigating pattern formation. Experiments in the 1980’s such as those conducted by 


Simonelli & Gollub (1989) tended to focus on the dynamics of the interaction of patterns 
in small containers excited by a sinusoidal forcing with a single frequency component. 
More recently the focus has switched to larger containers with multiple frequency com¬ 


ponents: as pointed out in Arbell & Fineberg (1998), one special feature of the Faraday 


experiment is that by using multiple frequency forcing one can investigate the interac¬ 
tion of a small number of controllable modes with different characteristic length scales. 
This provides some middle ground between pattern-forming systems where the initial 
instability involves just a single mode and in which complexity then results from sec¬ 
ondary bifurcations, such as in the Benard-Marangoni experiment or the Taylor-Couette 
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experiment, and cases where a system is highly turbulent and many spatial and temporal 
modes are excited. 

Within the Faraday experiment, a rich variety of complex patterns are seen, some of 
which have a complicated spatial structure but are time-periodic with the periodicity 
of the drive and some with both a complicated spatial and temporal structure. Aside 
from the ubiquitous stripes, squares and hexagons, observed patterns include: quasipat¬ 


terns (Christiansen et al. (1992), Edwards & Fauve (1994|); superlattice patterns (SLl) 
(Kudrolli et al. (1998); Epstein & Fineberg (2006)); spatially subharmonic superlattice 
states, modulated hexagonal disorder, two mode superlattices and unlocked states (Ar- 


bell & Fineberg (1998)); oscillons (Arbell & Fineberg (2000)), and double hexagon states 


|Arbell fc Fineberg (1998). An excellent summary of many of the experimental results is 


given in Arbell & Fineberg (2002). 


In parallel with the experiments, theoretical advances have resulted in an effective nu¬ 
merical method for performing the linear stability analysis of the Navier-Stokes equation 


that marks the transition from an unpatterned to patterned state (Kumar & Tuckerman 


(1994)). The challenge here is that many of the experiments are carried out at moderate 


viscosity whereas early theoretical results of Benjamin & Ursell (1954), which help ex¬ 


plain the underlying instability mechanism, are for an inviscid fluid. There have also been 
theoretical explanations for why, unlike the Benard-Marangoni system, the observed pat¬ 
terned states include superlattice patterns ([Silber fc Skeldon (1999); Silber et al. (2000); 


Silber & Proctor (1998); Rucklidge & Silber (2009)). 


These theoretical mechanisms rely on three-wave resonance between critical modes 
and modes close to critical and have successfully explained why particular superlattice 
patterns are observed. The suggested theoretical mechanism for the appearance of super¬ 
lattice patterns is compelling, and the link between which modes are excited and which 
patterns are observed has been explored in some detail experimentally, for example see 
[Arbell fc Fineberg (2002); Epstein & Fineberg (2006). However, without carrying out a 
careful quantitative comparison between experiment and theory, it is hard to know the 
extent to which the theoretical ideas really do explain the experimental results (Boden- 


schatz et al. (2000)). This is particularly true for the Faraday experiment where, as we 


will explain in greater detail below, the very region for which the theory predicts super¬ 
lattice patterns is the region for which some of the underlying assumptions of the theory 
break down. 

While it is not possible to write down closed-form solutions of the full Navier-Stokes 
equations for the Faraday problem, quantitative predictions of the patterns expected near 
the transition from non-patterned to pattern states can be made using weakly nonlinear 


analysis (Skeldon & Guidoboni (2007)). Weakly nonlinear analysis centres on using an 


asymptotic expansion in terms of the slowly varying amplitudes of the critical modes at 


onset, using ideas first developed in Malkus & Veronis (1958) and Segel & Stuart (1962) 


in the context of convection experiments. Analysis of the resulting amplitude equations 
leads to predictions on the relative stability of different patterned states. 

Our aim in this paper is to do a quantitative comparison between previously published 
experimental results and a weakly nonlinear analysis of the Navier-Stokes equations for 
the Faraday problem, with a particular focus on multiple frequency forcing. A previous 
comparison between the theory and experiments in the Faraday problem was carried out 


by Westra et al. { 2003) for single frequency forcing in which the authors declare that the 


excellent agreement they observe means that the Faraday problem is essentially a ‘solved 
problem’. In particular, their paper uses a Lyapunov stability argument based on the 


theory of Chen & Vihals (1999) to find the ‘most stable’ pattern as a function of the key 


non-dimensional groups in the problem representing non-dimensional measures of viscous, 
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gravitational and surface tension forces. Results from this theoretical study agree very 
well with the author’s experimental results. However, the extension to multiple frequency 
forcing is far from trivial. The method used by Chen & Vihals (19991 to reduce the Navier- 
Stokes equations to amplitude equations is not applicable, and the more general theory 


Skeldon & Guidoboni (20071 is needed; the addition of more frequency components 


introduces more non-dimensional parameters, resulting in a much greater ability to probe 
underlying three-wave mechanisms. Consequently, patterns such as superlattice patterns, 
which are not observed in the single frequency context, are found. In fact, as we will show, 
for multiple frequency excitation there remain many open questions. We note also that 


Westra et al. (2003) use a Lyapunov stability argument to determine preferred patterns. 


We use similar arguments, but have in addition carried out a bifurcation analysis of the 
relevant amplitude equations. This has the added benefit of not only determining the 
‘most stable’ pattern but also indicating regions where patterns are bistable. 

Specifically, in this paper we discuss to what extent the existing weakly nonlinear the¬ 
ory can explain observed patterns in multiple-frequency Faraday experiments and provide 
some new explanations in some cases. While agreement is very good in many cases, we 
note that quantitatively linking theory with experiment is particularly challenging after- 
the-fact as the results are sensitive to precise values of viscosity and surface tension and 
even the sign of the drive term, something that is not normally recorded. There are also 
places where the analysis strongly suggests that patterns should have a subharmonic 
component, when no subharmonic component has been observed. In particular, with re¬ 
gard to the superlattice patterns, we discuss two methods that have been used to promote 
the stability of superlattice patterns: firstly by approaching the so-called bicritical point 
in two-frequency forced experiments; and, secondly, by adding a third frequency to the 
drive. We highlight the differences between these two mechanisms. 

Overall, we confirm that three-wave interactions sit at the heart of why complex pat¬ 
terns are stablised. However, the discrepancies between theory and experiment suggest 
the need for further experimental and theoretical work to fully investigate the issues of 
pattern bistability and the role of bicritical/tricritical points in determining bifurcation 
structure. 


2. Equations 

Using a variety of container shapes, Edwards & Fauve (1994) elegantly demonstrated 
that many of the patterns that occur with moderate viscosity fluids in large containers 
are not strongly dependent on the lateral boundaries of the container. Consequently it 
is a reasonable modelling assumption to consider an infinite horizontal layer of viscous 
incompressible fluid of finite depth that is subjected to gravity g and to a vertical periodic 
excitation with non-dimensional frequency components ju). At the lower boundary the 
fluid is in contact with a rigid plane while at the upper boundary the surface is open to 
the atmosphere. This means that the upper surface is a free boundary whose shape and 
evolution is an unknown of the problem. 

The motion of the fluid can be described by the Navier-Stokes equations, where, to 
take account of the parametric excitation, a frame of reference which is moving with the 
periodic excitation is considered. The z-axis is chosen perpendicular to the rigid plane at 
the bottom, which lies at z = —h/l, where h/l is the non-dimensional depth of the layer 
when the fluid is at rest. A sketch of the geometry is shown in figure Assuming the 
free surface may be written as z = C(x,y;t), which excludes the formation of droplets 
or breaking waves, then the fluid motion in the bulk is described by the dimensionless 
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incompressible Navier-Stokes equations 

Vu = 0, 

dfU + u • Vu = — V7^ + CAu — (1 + f{t))e 3 , (2.1) 

where u = {u, v, w) is the velocity field, V the pressure and for multiple frequency 
excitation, 

m = E Gj cos{jujt + (2-2) 

j 

where j are integers and the non-dimensional amplitudes Oj and phases (pj are real. 

It is assumed that the bottom of the container is rigid so that at z = —h/l the fluid 
satisfies the no-slip boundary conditions 

u = V = w = 0. (2.3) 

At the free surface z = (^{x,y;t) we have the kinematic condition, which says that the 
surface is advected by the fluid, and two further conditions, one for the balance of the tan¬ 
gential stresses and one for the balance of normal stresses. This leads to three conditions 
at z = C(a:, y; t) namely 

dtC + udxC + vdyC = w, 

tx • Tn = t 2 • Tn = 0, (2.4) 

-r + 2C'nD(u)n = 3%- p^, 

where T = —VI + 2C'D(u) is the stress tensor, D(u) = (Vu -|- V^u)/2 is the rate-of- 
strain tensor, and % = Vh ■ (VijC/\/l + is the double mean curvature. Note 

that V = {'^H,dz) and V// = {d^, dy). The unit normal and tangent vectors are defined 
as 


n{x,y\t) 
ti{x,y\t) 
^2{x,y\t) 



The units of length, time, velocity and pressure have been taken as I, y l/g, gl and ggl 
respectively, where, g is the density of the fluid and g is the acceleration due to gravity. 
The length scale I is taken to be a length scale that is typical for the problem, such as the 
primary wavelength of the observed patterns. Here, Pe is the dimensionless pressure of the 
external ambient fluid and is assumed known. There are two non-dimensional parameters 
associated with the fluid, namely: C = vj{gPY/"^, the square of the inverse of the Galileo 
number, and B = d'/ggP, the inverse Bond number, where v is the kinematic viscosity 
and a is the surface tension. These two non-dimensional parameters measure the relative 
importance of viscous and surface tension forces compared to gravity respectively. There 
are three other sets of nondimensional parameters of importance, all associated with the 
excitation. These are: the non-dimensional amplitudes aj, frequencies jco and phases 4>j 
of the components of the excitation. In the comparison with experiments it is frequently 
useful to quote the relevant dimensional values: we have used the convention that all 
dimensional variables are labelled with a tilde. 
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Figure 1. Sketch of a cross-section through the layer of fluid. 


As in Kumar & Tuckerman (1994), it is convenient to define a new pressure, 

p = V + il + fit))z, 


(2.5) 


which has the effect of shifting the acceleration term from the momentum equation to 
the normal stress condition. In addition, we eliminate the pressure from the momentum 
equation by taking —(V x Vx). Using the relation VxVxu = V(V-u) — Au and the 
fact that V • u = 0, the problem then becomes 


V • u = 0, 

dtAu — CAAu = V X V X (u • Vu), 
with boundary conditions on z = —h/l, 

u = V = w = 0 


and on z = C, 


( 2 . 6 ) 

(2.7) 


dtC + udxC + vdyC = w, 

ti • Tn = t 2 • Tn = 0, 

2CnD(u)n = BH + p — Pf. — {I + 


( 2 . 8 ) 


Equations (2.6) with boundary conditions (2.7) and (2.8) have a trivial solution, 

u = 0, p = pe, C = 0- 


(2.9) 


This solution corresponds to an unpatterned state where there is no relative motion of 
the fluid with respect to the moving frame so the surface of the fluid is flat. 


3. Does linear theory agree with experiments? 

Understanding the linear stability of the unpatterned state is at the heart of under¬ 
standing many of the nonlinear patterns that are formed close to onset. This is because 
when the unpatterned state becomes unstable to one critical mode with a given critical 
wavenumber, there are often several other modes with different wavenumbers that are 
themselves only weakly damped. This is particularly true with multi-frequency forcing. 
Resonant interaction of instabilities from different critical/close to critical modes drive 
the selection mechanisms for the occurrence of particular patterns (Silber et al. (2000)). 
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Benjamin & Ursell (1954) recognised that the linear stability of the flat-surface solution 


for an inviscid, infinite depth fluid driven by a single frequency reduces to a Mathieu 
equation. The Mathieu equation contains two parameters, related to the frequency w and 
amplitude d of the excitation respectively. Solutions to the Mathieu equation divide the 
parameter plane into regions of bounded and regions of unbounded solutions where the 
regions of unbounded solutions form tongues that touch the frequency axis at frequencies 
muj/2,m = 1,2..., the largest tongue occuring for m = 1. The tongues are typically 
classified as either harmonic or subharmonic, depending on whether or not they are 
an integer multiple of the frequency w. This picture is modified with the addition of 
damping: the boundaries of the regions are perturbed and no longer touch the frequency 
axis, consequently a finite amplitude of excitation is required to excite waves; the regions 
of bounded solutions become regions where the unpatterned state is locally stable; the 
unbounded regions become regions where the unpatterned state is unstable. 


Kumar & Tuckerman (1994) identified a numerical method to find the instability 


tongues that can be used for all fluid viscosities and all depths and applied it to the 
case of single frequency excitation. This was extended to multiple frequency excitation 
by Besson et al.\ (|1996[). A typical example of tongues computed using Besson et al. 


(1996)’s method and the corresponding bifurcation set for the primary stability bound¬ 


ary is shown in figure 

The linear stability problem appears to be solved: the numerical method works well and 
Besson et aL\ (|1996[) the authors show that there is excellent agreement between their 


linear stability calculations and experiments. However, unless experiments and theory 
are carried out hand-in-hand, in practice there remain some difficulties in obtaining 
really good agreement between theory and experiment even for this first transition from 
patterned to unpatterned state. This is illustrated in figure where numerical linear 
stability results are superimposed on the experimental bifurcation sets published in: 

( |1994 ) (panels (a) and (b)); Kudrolli et al.\ (|1998 ) (panel (c)); Epstein 


Edwards & Fauve 


& Fineberg (2006) (panel (d)), and Ding & Umbanhowar (2006) (panels (e) and (f)). 


The central issue here is that the position of the curves is sensitive to the values 
of the surface tension, density and the viscosity and yet the values quoted in papers 
are often taken from the manufacturers’ specifications for the fluids used. In order to 
illustrate the issue, in each case we have plotted the linear stability curves for the quoted 
viscosity and the quoted viscosity plus or minus 5%, a typical quoted tolerance for the 



the actual viscosity of the fluid was higher than the upper value (panel (d)). We have 
focussed here on viscosity because for the different ffuids used in the results presented 
here, errors in the viscosity have the biggest effect on the linear stability boundary. An 
error in the viscosity of 5% can result in an error in the linear stability boundary of 5%, 
whereas a 5% error in either the surface tension or the density leading to a 5% error in 
the ratio a/g leads to an error in the linear stability boundary of only around 1%. 

Our aim in this paper is to compare weakly nonlinear theory with experiments. Since 
the transitions from one preferred patterns to another turn out to be quite delicate, 
clearly it is only realistic to hope for agreement for nonlinear pattern selection if there is 
first excellent agreement with linear theory. In making our nonlinear comparisons below, 
we have therefore used values of the viscosity found by fitting the linear theory to the 
published data. 
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Figure 2. Ty pical tongues and lin ear stability region for the unpatterned st ate computed using 
the method of Besson et al. (19961. Parameters as for Kudrolli et al. (19981: the (dimensional) 
excitation is f{t) = aecosGoit + fiy cos(7a)t + (f)?), where dg — acosx, 3? = asiny, (pj = 20°, 


wj^-K = 16.44 Hz, a = 20.6 dyn cm 


Q = 0.95 g cm °, h = 0.3 cm. Note 


that the structure inherited from the Mathieu equation is seen by the fact that the minimum 
of each tongue occurs at half integer multiples of the drive frequency. In (a) x — 0°, the forcing 
reduces to /(t) = ag cos 6ujt, with a forcing frequency of OlO. The largest tongue has a minimum 
for wavenumber of approx 15cm“'^ for a forcing amplitude d of approx 7.5g and corresponds 
to a mode of frequency In (c), y = 90° the forcing reduces to f{i) = drcosltoi with a 
forcing frequency of 7Cj leading to a primary instability mode with frequency 7ci)/2. In the case 
of (b), X = 63°, the presence of both 6a) and 7a) components in the forcing frequency mean 
that the drive has periodicity 27r/a) and is therefore of frequency a). This leads to tongues at 
ma)/2 ... ,m = 1,2,.... The first five tongues that are visible from left to right correspond to: 
a)/2,4a)/2, 5a)/2, 6a)/2, 7a)/2. The largest tongues correspond to the 6a)/2 and the 7a)/2 modes 
and are driven by the two main frequency components of the drive. There are tongues corre¬ 
sponding to modes with frequency 2ub/2 and 3a)/2 but these occur off the top of the region 
shown. This particular value of x is close to the ‘bl-critical’ point where both harmonic modes 
with frequency 6a)/2 and sub-harmonic modes with frequency 7a)/2 onset simultaneously. It is 
close to this bi-critical point that many of the exotic patterns are observed, (d) Bifurcation set 
showing the position of the tongue minimum that marks the instability of the unpatterned state, 
as a function of fig and dr. 


4. Predicting patterns close to onset 

4.1. Theoretical ideas 

Once a mode with a given wave number has become unstable, the fact that there is no 
preferred horizontal direction means that standing waves of any orientation can occur. 
In the Benard-Marangoni experiment, which has a similar orientational invariance, this 
means that patterns consisting of nonlinear superpositions of modes with wave vectors of 
















Figure 3. Comparison of the experimentally measured transition from unpatterned to pat- 
terned state with numerically calculated linear stability curves using the method of |Besson| 


et al. (19961. For each of the numerical calculations, three different values for the vi scosity are 
used, (a) arid (b) Nu merical curves overlaid on figure 12 of Edwards & Fauve (|1994 1 . Adapted 
„,ni, Edwards & Fauve _ cc „™-l. _ i no _l a AC „™2 „-i. 


with permission from 
,-3 


(19941. d = 65 dyn cm u = 1.00 ± 0.05 cm” s 
g = 1.22 g m“'^; h = 0.29 cm; j = {4,5}; 4>a = 0; 0 5 = 75°; = 14.6 Hz. (c) Numerical 

curves overlaid on figure 6(a) of Kudrolli et al. (19981 reprinted from Physica D, Vol 123, A. 
Kudrolli, B. Pier and J.P. Gollub, ‘Superlattice patterns in surface waves’, p99-lll. Copyright 
(1998), with permission from Elsevier, a = 20.6 dyn cm“^; z> = 0.20 ± 0.01 cm^ s“^; g = 0.95 g 
m“®; h = 0.3 cm; a,-, ) = |4, 5|; 04 = 0; (j)^ = 16°; — 22 Hz. (d) Numerical curves overlaid 


on figu re 2 of Epstein & Fineberg (20061. Adapted with permission from|Epstein & Fineberg 
(20061. Copyrighted by the American Physical Society, a = 20.6 dyn cm“ , = 0.180 ± 0.009 

cm s“^; g = 0.949 g m“®; h = 0.3 cm; aj,j = |6,7}; <j)i = 0°; uiI'I'k = 14 Hz. (e) and (f) 
Numerical curve s overlaid on figures 1 and 3( a) of Ding & Umbanhowar (20061. Adapted with 
permission from [Ding fc Umbanhowar| ( |2006| ). Copyrighted by the American Physical Society, 
d = 20.6 dyn cm“ ; v = 0.20 ± 0.01 cm s“^; g = 0.95 g m“®; h — 0.65 cm; aj,j = {4,5}; 

= 0; <))5 = 16°; Cb/2'K = 20 Hz. 
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Figure 4. Possible spatial three-wave resonant interactions defined such that ki -I- k2 = ks. 
(a) |ki| = Ikal = Ikal; (b) |ki| = jksl > [kal; (c) |ki| = jkal < IksI; (d) |ki| / jkal / |k3|. 


the same wave number but different orientations often lead to the observation of patterns 
with, for example, hexagonal symmetry. In principle, superlattice patterns could also 


occur in convection experiments (Skeldon & Silber (19981), but in practice they have not 


been seen without the addition of a vertical oscillation, as in Rogers et al. (2000). 


The ubiquitous occurrence of hexagons is a consequence of the importance of three- 
wave interactions in determining which patterned states occur: three-wave interactions 
give rise to the lowest order nonlinear (quadratic) terms in the amplitude equations that 
describe behaviour close to onset. For Benard-Marangoni convection, where instability to 
a single dominant wavenumber occurs, the three-wave interaction of importance occurs 
between three wave vectors with the same critical wavenumber, as shown in figure Qa). 
A distinct feature of the Faraday problem is that the Mathieu tongue-like structure to 
the linear stability problem means that although typically there is a single mode that 
becomes unstable first, there are nearby modes with different wavenumbers that are only 
weakly damped. This can give rise to other three-wave interactions of relevance, such as 
those shown in figure |^b) , (c) and (d). 

These three-wave interactions involve waves with wavevectors ki, k 2 and ka and re¬ 
spective frequencies nuj(2,pui j2 and quj. For the three waves to interact they must satisfy 
a spatial resonance condition, 

±ki±k2=k3, (4.1) 

and a temporal resonance condition. In the spatial resonance condition the choice of 
sign arises because the waves are standing waves and so have spatial fourier components 
with both signs of wave vector. The temporal resonance condition depends on both the 
frequencies of the waves and the various frequencies contained within the forcing term 
f{t) and requires 


n p 

* 2*1 


±j = q 


(4.2) 


where j is one of the frequency components of the drive see equation (2.2 ),cf Topaz et al. 


(2004) and where all possible sign combinations of the different terms on the left-hand 


side need to be considered because the waves are standing waves and the forcing is real. 

The fact that a temporal resonance condition needs to be satisfied is one feature that 
distinguishes the Faraday problem from Swift-Hohenberg multiple resonance problems 
such as those studied in |Miiller (1994), Lifshitz & Petrich (1997), Rucklidge et al. (2012) 
and others. This means that results from Swift-Hohenberg-like equations need to be 
interpreted with care when applied to the Faraday problem, a point that we will return 
to in the discussion. 

For the appearance of superlattice patterns the argument goes that, given three critical 
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modes kj as shown in figure Qb) or (c) with amplitudes Ai, then using a multiple 
timescale expansion near onset would lead to equations for the evolution of the amplitudes 
on a slow timescale of the form 


Ai = Ai^i + 01^2^3 + j4i(a|Aip + 6o|Al2p + + ... 

Ai = A 1 AI 2 + 01^11^3 + A2(&oI^iP + + c|A 13 |^) + . . . 

A 3 = A2^3 + 02^1^2 + ^3('^|^l|^ + d|^2|^ + + ■ • ■ , 


(4.3) 


where Ai and A 2 are the linear growth rates of the respective modes, and oi, 02 , a, bo, c, d 
and e are all real-valued constants. The quadratic coefficients, ai and are non-zero 
only if the temporal resonance condition is met. Note that t here is no assumption of 
weak forcing/damping and so terms such as Ai do not appear (Alnahdi et ai (2014)). 

Now, the mode with wavevector ka is not at its critical point but is weakly damped, 
so A 2 < 0 and, close to onset (i.e. |Ai| small enough), it is therefore slaved by the critical 


modes. Consequently, one can perform a centre manifold reduction on equations (4.3), 
resulting in 


and 


where 


^3 = -^^1242 
A 2 


Ai — XiAi + Ai{a\Ai\'^ -\- &(0res)l^2p) + • ■ ■ 

A 2 = A 1 A 2 -|- A2(&(0res)|4llp -|- a|A2p) -)-..., 


^(^res) — ^0 T 


bres — 


Q!iQ!2 

Aq 


(4.4) 


(4.5) 


The presence of the weakly damped mode therefore changes the value of the cross¬ 
coupling coefficient 5(0res) between modes with wavevectors ki and k 2 , offset at an 
angle 0res- The value of 0res is determined by the ratio of critical and weakly damped 
wavenumbers, but the same idea holds for any three-wave interaction between any ki 
and k 2 on the critical circle. This results in a function b(0) with either a distinctive peak 
or dip at 0 = depending on the sign of aia 2 and on the weakness of the damping 
for the weakly damped mode. 


Analysing the amplitude equations (4.4) shows that there are two types of solutions 


that bifurcate from the trivial solution Ai = A 2 = 0, namely stripes (Ai ^ 0, A 2 = 0 
or vice versa) and rectangles (Ai = A 2 ). The relative stability of rectangles to stripe 
perturbations is dependent on the relative size of the self-coupling coefficient a in equation 
(4.4) and &(0res) where, if bres > 0 (aia 2 > 0), then the stability of rectangles is enhanced 
by the three-wave interaction and if bres < 0 (0102 < 0), then the stability of rectangles 
is suppressed. 

Of course, in the Faraday problem there are not just two modes and there are many 
interactions and several weakly damped circles, but this idea that three-wave resonances 
can promote patterns associated with the angle 0res is powerful (Silber et ai (2000)). The 
idea is that the dispersion relation determines which wavenumbers are critical or close to 
critical; allowed three-wave resonant interactions then select out particular wavevectors; 
the allowed wavevectors are necessarily oriented at particular angles as determined by the 
ratios of available wavenumbers, leading to specific values for Ores- Superlattice patterns 
essentially consist of a nonlinear superposition of two sets of hexagons offset at an angle 


to each other. There is a whole family of different superlattice patterns (Dionne et al. 


(1997)), each corresponding to a different angle, but the particular superlattice that will 


be promoted will be that with an angle given by 9r 
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The superlattice patterns observed by Kudrolli et al. (1998) for {6,7} excitation ex¬ 
emplify the idea and motivated the early theoretical work. The superlattice patterns 
that they observe are spatially periodic and a spatial Fourier transform of the pattern 
indicates that they consist of essentially two sets of hexagonal modes offset by an angle 
of 22°. In this case, the particular three-wave resonance of relevance is of the type illus¬ 
trated in figure l^b) and is between a harmonic mode that is weakly damped of frequency 
2 J)/ 2 , related to the 2 a )/2 tongue, and two wavevectors with wavenumbers corresponding 
to the 6(h/2 tongue ( [Silber et al. ( 2000| ). This was surprising because the superlattice 
state appears near the bicritical point where both QQ/2 and 7Q/2 onset simultaneously 
and naively, from inspection of the linear stability diagram (see figure [^b)) one would 
assume that it was a result of the interaction of 6 w /2 and 7ui/2 modes—the 2uj/2 tongue 
onsets at an amplitude of approximately 50 so is not even visible on the scale shown. 

This argument not only explains the presence of the 22° superlattice patterns but also 
why patterns with angles close to 30° are seen in j = {4, 5} forcing (Skeldon & Guidoboni 


(2007)). Extensions of this basic idea have been used to suggest ways to design forcing 
frequencies to promote particular patterns in both Faraday waves (Porter et al. ( 2004| )) 
and in a model partial differential equation (Rucklidge & Silber (2009)). 


The results of Kudrolli et al. (1998) and the theoretical results of Silber et al. (2000) 


showed how modes visible in the spatial Fourier transform of the experimentally ob¬ 
served pattern can be linked to critical/weakly damped modes and three-wave resonant 
interactions. It follows that given an experimentally observed pattern it should be pos¬ 
sible to confirm whether or not it is a result of a particular interaction by establishing 
which modes are present in the pattern. However, identifying precisely which modes are 
involved can be tricky, as we illustrate in the following section. 


4.2. Practical identification of relevent modes 


In Epstein & Fineberg (2006) a detailed comparison of different patterns, referred to as 
‘grid’ states and labelled as 3 : 2, 4 : 3 and 5 : 3 respectively is presented. The spatial 
Fourier transform of the patterns enables the identification of the wavenumbers that are 
present in these patterns: in each case they appear to be dominated by modes sitting on 
two different circles, one on the primary harmonic instability with wavenumber kc and a 
second circle with respective wavenumbers ^ 3,2 = 0.38fcc) ^ 4:3 = O.SSfcc and ^ 5,3 = 0.23fcc. 
These ratios are fixed by the fact that the observed grid states consist of modes that fit 
on to a hexagonal lattice. This suggests that they are strong candidates for superlattice 
patterns that occur as a result of the mechanism discussed in section 4.1. |Epstein~^ 
Fineberg (2006) computed the linear stability curves, and their linear stability diagrams 
are re-computed in figure]^ a) and (b ). We note that the linear st ability curves are similar, 
but not identical, to those given in Epstein & Fineberg (2006) particularly at small k. 


This is because we have made a different choice in the sign of the forcing term to that 
made by Epstein^ fc Fineberg (2006). This is explained further in Section]^ As stated by 
Epstein & Fineberg (2006) and, as can be seen in this figure, the ^ 3,2 wave number lines 


up with the first harmonic tongue (2Cj/2) and this strongly suggests that it is the same 
i necha nism that produces the familiar 22° superlattice patterns, seen by [Kudrolli et aL\ 
{ 1998 1), that produces the fc 3.2 g rid state. For the two other states the situation is less 


clear. [Epstein fc Fineberg (2006) state that ^ 4,3 lines up with the the second harmonic 
ton gue {Aw 12'). but this is not c lear from figure]^ (the position of the line drawn in figure 


5 of 


Epstein fc Fineberg (2006) is not consistent with the value of 0.55fcc given). Epstein 


fc Fineberg (2006) observe that the ^ 5,3 mode does not appear to be aligned with any 


tongue and that the 5 : 3 pattern occurs in a wedge-shaped region of parameter space 
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(a) 


(b) 





Figure 5. Linear stability curves and Floquet multipliers calculated from the finite depth 
Navier-Stokes equations fa~) and fc) h — 0.3 cm (b) and (d) h = 0.2 cm. Parameter values 
are taken from 


Epstein & Fineberg (|2006|) and are a = 20.6 dyn cm v = 0.18 cm^ s 

y = 0°; £u/27r = 14Hz. 


g = 0.949 g m dj,j = { 6 , 7, 2}, where de = 3.4 ,07 = 6.4, 0,2 = 0.5; 


that emerges from the bicritical point, unlike the 3 : 2 and the 4 : 3 state that occur at 
a short distance from this point. 

The difficulty here is that the linear stability diagram indicates that at an amplitude 
of excitation of approximately 7.5g the mode with wavenumber kc and frequency 6a;/2 
becomes unstable and that there are five other damped modes with smaller wavenumber. 
But the position of the tongues for these modes only gives a rough idea of the precise 
wavenumber and damping associated with each mode at the pattern onset amplitude of 

7.55. 

More accurate information on the damped modes can be obtained by calculating the 
most critical Floquet multipliers, as detailed in Appendixj^ The results for the parameter 
values used by Epstein & Fineberg (2006) are shown in figure [^c) and (d). We note that 
with either p ossible sign for the forcing, the Floquet multipliers look very similar. As 
identified by Epstein & Fineberg (2006), the ^ 3,2 mode is aligned with the harmonic 


tongue associated with frequency 2uj/2 for h = 0.3 cm, but not for h = 0.2 cm and 
supports their observation that the k^-^ grid state will be seen for h — 0.3 cm and not 
for h = 0.2 cm. The remaining two wavenumbers ^ 5,3 and ^ 4,3 are very close to the 
subharmonic wa ves of frequencies 6j/2 and 3 uj /2 respectively. This clarifies the issue that 
the ^ 5,3 mode in Epstein & Fineberg (2006) did not appear to align with any mode on 


the linear stability graphs and strongly suggests that the fc 4:3 mode was mis-identified as 
harmonic. However, it also presents a problem in that temporal constraints mean that 
it is not possible for two harmonic wavevectors to form a three-wave resonance with a 


single subharmonic wavevector. Consequently, ai = 02 = 0 in equations (4.3), so the 


discussion in section 4.1 does not hold as it stands and there is no special angle 0^ 
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Epstein & Fineberg (2006) point out that the ^ 5:3 case could be a result of a three- 


wave resonance between modes with frequencies Co/2, 6w/2 and 7Co/2 and this is strongly 
supported by our Floquet multiplier calculation and is a permitted resonance in that it 
satisfies both spatial and temporal constraints. Such a resonance appears to be consistent 
with the spatial Fourier transform in figure 3(c) of Epstein & Fineberg (2006) and is 


similar in structure to previous states identified in figure 20(c) of Arbell & Fineberg 


(2002) as 2MS states. This would make it a three-wave resonance of the form shown 
in figure (d). This kind of resonance was discussed briefly by [Porter fc Silber (2002) 
and results in patterns that occur in a wedge that emerges from the bicritical point. 
This is consistent with the bifurcation sets shown by Epstein and Fineberg. The only 
remaining conundrum for this particular pattern is that Epstein and Fineberg state that, 
although all the evidence points to subharmonic modes being important, no subharmonic 
component was found. 


5. Observed patterns near onset: theoretical bifurcation sets 
compared with experiment 

In spite of the difficulties in identifying modes, nevertheless, in experiments with dom¬ 
inant forcing frequencies j = {4,5} or j = {6,7} there is a coherent picture: in the 
{ 6 , 7} case, observed superlattice patterns near onset have an angle of 22°; for j = {4,5}, 
observed patterns near onset are quasipatterns with an angle of 30°. Although there 
are technical difficulties with considering quasipatterns in the same way as superlattice 


patterns (Rucklidge & Rucklidge (2003); Rucklidge & Silber (2009); looss & Rucklidge 


( 2010 )), nevertheless they appear to fit within the same framework with resonant in¬ 
teractions with weakly damped modes associated with 2Co/2 contributing to bres and 
explaining the appearance of the appropriate angle. The different angle seen in the {4, 5} 
case as compared with the { 6 , 7} case is because of differences in the wavenumbers of 
the interacting modes. In both cases, super lattice/quasipatterns are only seen near the 
bicritical point where the two modes driven by the two main components of the forcing 
onset simultaneously. 

Here, we aim to carry out a careful comparison of the results of weakly nonlinear the¬ 
ory and experimental results to explore to what extent there is quantitative agreement 
between experiments and theory as the bicritical point is approached. We present results 
from the experiments of Ding & Umbanhowar (2006) because we have excellent agree¬ 
ment for the linear stability curves with v = 0.21 cm^ s“^, as shown in figure [^e) and 
(f), and this set of experiments includes the most comprehensive study of how multiple 
frequencies interact via the superlattice mechanism outlined above. Specifically, acknowl¬ 
edging the importance of the 2Co/2 tongue. Ding & Umbanhowar (2006) carry out a series 
of experiments that systematically explore the effect of adding a third frequency that pro¬ 
motes this mode. 

The weakly nonlinear theory is carried out using the method given in [Skeldon fc| 


Guidoboni (2007) for analysing equations (2.6) to (2.8). For all details of how weakly 


nonlinear analysis is used to derive coefficients for the relevant amplitude equations we 
refer the reader to Skeldon & Guidoboni (2007). In order to examine the relative stability 
of superlattice patterns a minimal set of twelve amplitude equations corresponding to 
twelve vectors on the critical circle are needed. However, the stability of the planforms 
that bifurcate from the non-patterned state can be found by calculating the coefficients 
on three one-dimensional subspaces of this twelve dimensional problem: one describing 
stripes, one for hexagons and one for rectangles generated from wavevectors separated 
by an angle 9 as detailed inlSkeldon & Silberl (|1998p. By varying 9, this formulation then 
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Figure 

6 

7 

5 

g 

TO 

TT 

T2 


Table 


{4,5} (04,</>5) = (O,16°) 

{4, 5, 2} fla = 0.85, (04, 05. 02) = (0°, 16°, 32°) 
{4,5,2} 55=53, (04,05,02) = (0°, 5°, 32°) 

{4, 5, 2} (54,55, 52) = (3.8, 5.8, 0.8)3, <(>4 = 0° 

{6,7 (06,07) = (0°, 40°) 

{6, 7, 2} 57 = 73, (06, 07. 02) = (0°, 40°, 80°) 

{6, 7, 2} (56, 57, 52) = (5.2, 7.8, 0.6)3, < f>e = 0° 


Parameters 

varied 

54 and 55 
04 and 05 
04 and 0.2 
05 and 02 
06 and 3.7 
oq and 02 
07 and 02 


Figure from 


Ding & Umbanhowar 

- 1 - 


(20061 


3 (b) 

4 (a) 
4(b) 

2 

6(a) 

6(b) 


1. Summary information for f igures comparing weak l y non linear analysis with the 
experiments of Ding & Umbanhowar (20061. 


results in determining the stability of stripes to any perturbation of the same wavenum¬ 
ber; the stability of hexagons to (i) stripes arbitrarily close to any given orientation (ii), 
rectangles arbitrarily close to any given aspect ratio, (iii) to superlattice patterns of any 
angle; the stability of rectangles to (i) stripes of any orientation (ii) superlattice pat¬ 
terns on the same lattice; the stability of superlattice patterns to hexagons, stripes or 
rectangles that are made up of a subset of the superlattice wavevectors. 

Since the weakly nonlinear stability analysis indicates that for many parameter values 
there is more than one stable pattern, results are presented in two ways, first showing 
the regions of stability for individual patterns—only those that have some region of sta¬ 
bility are shown, and secondly showing the bifurcation set for the most stable state, as 
computed from the weakly nonlinear coefficients using the Lyapunov functional given in 


Skeldon & Guidoboni ( 

2007 

). The most stable states are superimposed on the experi- 

mental results of 

Ding & Umbanhowar 

(2006 

). In the case of superlattice patterns, the 


resonant interactions mean that it is always the particular pattern associated with Ores 
that has the largest region of stability and it is perturbations associated with this angle 
that first destabilise hexagons. Consequently where superlattice patterns are shown it is 
always the superlattice pattern associated with Ores that is relevant. 

Ding and Umbanhowar consider the two sets of frequency components most widely used 
by other experimental groups, namely j = {4, 5} and j = { 6 , 7} and then systematically 
investigate the inclusion of a third mode, so considering the combinations j = {4, 5, 2} 
and j = { 6 , 7, 2}. We have carried out the weakly nonlinear stability analysis for each of 
their parameter studies and the results are presented in figures 6 p 2 where the principle 
parameters considered and the corresponding figures in Ding & Umbanhowar (2006) are 
summarised in table [T] 

In each case the weakly nonlinear results are superimposed on the experimental re¬ 
sults. The parameters have been chosen to match those quoted in the paper of |Ding ^ 
Umbanhowar (2006) with the exception of the viscosity where, as seen in figure the 
linear stability curves fitted best with a viscosity ofn = 0.21 cm^s“^ rather than z> = 0.2 0 
cm^s“^ quoted in their paper (or the value oi i> = 0.204 cm^s“^ quoted in Ding (2006)). 


5.1. Results for {4,5} and {4,5,2} 

In figure we see that the weakly nonlinear theory predicts bistability between squares 
and hexagons for low values of 05 , with squares losing stability as the bicritical point 
is approached. There is good agreement between the point at which squares become 
unstable (approximately as = 2.2g, figure [^a)), and the transition between squares 
and hexagons in the experiment. However, before squares become unstable, the theory 
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05(5) 


Figure 6. Comparison with figure 1 from 
di = {4,5},(<?i4,</>5) = (0,16°),(h/27r = 20 Hz. The 


Ding 

fluid 


& Umbanhowar (20061. 


tion to figure 3(e),(f) with v — 0.21 cm'^s . (a)-(b) 


parameters are as m the cap- 
Stability regions for different planforms. 


(a) Stable squares in green; unstable squares in grey; (b) stable hexagons in blue; unstable 
hexagons in grey. Note that for as > 0, there are quadratic terms in the amplitude equations 
which mean that hexagons bifurcate transcritically. The unstable solutions that are produced 
in this transcritical bifurcation are stabilised in a saddle-node bifurcation, resulting in a 
small region of stable, subcritical hexagons. This region of stable subcritical hexagons is only 
visible on the scale of the figure close to the bicritical point, (c) ‘Most stable’ state using a 
Lyapunov energy argument: -f squares; * hexagons; □ q uasipatterns. The theoretical r esults 
are overlaid on figure 1 from Ding & Umbanhowar (20061. Adapted with permission from Ding 


& Umbanhowar ( 2006[). Copyrighted by the American Physical Society. 


suggests that there is a region of bistability between hexagons and squares and that in 
this region hexagons and not squares, as seen in the experiment, are the most stable state. 
Of course, in this region, the particular pattern observed experimentally will depend to 
some extent on the way in which the experiments were carried out. For example, if the 
experimental procedure was to fix 04 and to increase 05 from zero in small steps, one 
would expect to see squares until the limit of their stability. Regions of quasipatterns are 
also predicted by the weakly nonlinear theory (almost off the scale of the figure and 
observed in the experiment but since they do not occur close to onset it is perhaps not 
surprising that the theory and experiment do not agree. 
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(b) 


3 



asig) 



Figure 7. Comparison with 
j = {4,5,2},(04,^5,</>2) = (0,16' 


Ding & Umbanhowar (20061) Figure 3(b), 
32“). (a)-(b) Stabiiity regions for different planforms. 
(a) Stable hexagons in blue, unstable hexagons in grey; (b) stable quasipatterns in red, 
unstable quasipatterns in grey. Both quasipatterns and hexagons bifurcate transcritically 
and may be stabilised in a saddle-node bifurcation. This can result in stable subcritical 
quasipatterns/hexagons, (c) ‘M ost stable’ planform using a Lyapunov energy argument 
superimposed on figur e 3(b) from Ding & Umbanhowar (20061. Adapted with permission from 
[Ding &: Umbanhowar (20061. Copyrighted by the American Physical Society. * hexagons; □ 
quasipatterns, • disordered states, where the coloured symbols are the theoretical results and 
the black symbols the experimental results. The red line is the linear instability curve for the 
harmonic 2a}/2 tongue. The black cross indicates the point at which the experiment shown in 
figure [^c) was carried out. 


In fig ure [7| we see that by a dding in a forcing component that excites the 2Cj/2 mode 
directly Ding & Umbanhowar (2006) showed that they could increase the region of quasi¬ 


patterns, to the extent that they become the first pattern observed after the unpatterned 
state becomes unstable. Our theoretical results for the same parameter values also show 
an enhanced region for the superlattice patterns—but not as enhanced as for the exper¬ 
iment. Two possible reasons for this are (i) the nearby presence of the bicritical point 
where both the 2uj/2 and the 4d;/2 tongue onset simultaneously which means that our 
weakly nonlinear calculations would need to be extended; (ii) the sensitivity of the results 
to phase. We discuss both of these in more detail below. 

Figurej^shows excellent agreement between the theoretical results and the experimen¬ 
tal results for the pattern at onset for increasing 0 , 2 , with agreement diminishing with 
distance from onset. 

There is surprisingly good agreement for the variation with phase shown in figure 
In the experimental results, black indicates low correlation at an angle of 30°. This 
lines up well with areas where the theory indicates competition between squares and 
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02(5) 


0,2 (q) 


(C) 



Figure 8. Comparison with |Ding & Umbanhowar| (|2006|) figure 4(a). 
j = {4, 5, 2}, 05 = 5g, {(j)A, 4>5, <(> 2 ) = (0,16“,32“). (a) and (b) stability regions for hexagons and 
quasipatterns respectively. Coloured symbols stable; grey unstable. Hexagons and quasipatterns 
both bifurcate from the unpatterned state at a transcritical bifurcation, where the unstable 
bifurcating branch changes direction at a saddle-node bifurcation. This leads to regions of 
stable hexagons and quasipatterns that are subc ritical. (c) ‘Most stable ’ planfo rm: * hexagons; 
□ quasipatterns, overlaid on figure 4(a " ~ 


(a) from Ding & Umbanhowar| ( 2006|). Adapted with 
( |2006 1 . Copyrighted by the American Physical Society. 


permission from Ding & Umbanhowar 
The coloured symbols are the theoretical results and the black symbols the experimental 
results. The black dots indicate disordered states. Note that the linear stability boundary for 
the 2u3/2 mode is just off this diagram, and is approximately parallel to the a 4 axis through a 
value of a 2 ~ 1. The bicritical point where the 2ih/2 and the 4u}/2 modes onset simultaneously 
occurs when ( 22 , 04 ) = (1.05, 3.80)g. 


hexagons. White areas, where there is high correlation at 30°, line up well with stable 
quasipatterns. In the regions where there are no coloured symbols, the weakly nonlinear 
analysis indicates no stable states. This region is mostly mid-grey in the experiment, 
indicating some correlation at 30°. 
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The agreement in figure is surprising because the results correspond to a point with 
04 = 3.8g, 05 = 5.9g and 02 = 0.8g, which according to the theory is close to a transition 
from quasipatterns to hexagons in a region where both hexagons and quasipatterns are 
stable, whereas in the experiment it is firmly in the quasipattern region. (In order to help 
cross-reference between the two figures, a black cross has been marked on both figurej^c) 
and figure]^) Figure [^highlights the sensitivity of the results to phase and in fact, with 
our weakly nonlinear analysis we find that the size of the region of quasipatterns changes 
significantly depending on the phase. This sensitivity to the precise value of the phase 
could be part of the difficulty in obtaining agreement for the onset of quasipatterns in 
figure 

The results are, of course, also sensitive to the viscosity. Away from the bicritical point, 
changing the viscosity tends to translate the nonlinear bifurcation lines as one might 
expect: upwards if the viscosity is increased, downwards if the viscosity is decreased, in 
line with the changes to the transition from unpatterned to patterned state seen in Figure 
3. Close to the bi-critical point the effects are more pronounced with small changes in 
viscosity sometimes leading to large changes in the regions in which particular patterns 
dominate. 


5.2. Results for {6,7} and {6,7,2} 

In figurej^the weakly nonlinear theory predicts bistability between squares and hexagons 
for low values of dy, with squares losing stability as the bicritical point is approached. As 
for figure although there is good agreement between the point at which squares become 
unstable (approximately 07 = 4.5g), and the transition between squares and hexagons in 
the experiment, this agrees less well with the results of the energy argument. 

In figure both theory and experiment agree that for low and moderate values of 02 
hexagons are the preferred pattern. As 0,2 increases, hexagons are replaced by superlattice 
patterns, although the transition observed experimentally occurs slightly earlier than in 
the weakly nonlinear calculations. Away from onset, the agreement is qualitative rather 
than quantitative: in the experiments superlattice patterns are observed for a much larger 
region than predicted by the weakly nonlinear theory and hexagons are observed in 
regions where the theory predicts bistability of rectangles and hexagons, with rectangles 
being the most stable state. For large values of 02 and de, the experiments see disordered 
states. The weakly nonlinear theory cannot predict such states, but note that none of 
the planforms considered are found to be stable for large 62 and dg (there are no symbols 
from the theoretical calculations in the top right of figure |TT|) . 

The results shown in figure again show that the theoretical results have the same 
diagonal dependency as seen in the experiments, a consequence of a phase invariant as 
discussed further below. For the experimental results, white indicates a high correlation 
with the angle of 22 ° and this does seem to align with where superlattice patterns occur, 
although for much of the lighter regions the weakly nonlinear analysis predicts no stable 
pattern. The black regions correspond to a low correlation at 22° and |Ding fc Umban-| 
howar (2006) state that this region contains ‘disordered’ patterns: our theoretical results 


suggest that both rectangles and hexagons are stable for much of this region each having 
a similar Lyapunov energy. Consequently one might expect competition between these 
two states, resulting in the observed disorder. 


5.3. Common themes 

From our comparison, a number of themes emerge: 

( а) There is excellent agreement with theory at onset in some cases (figure . 

( б ) In other cases, figure]^ and 10 bistability of patterns makes it difficult to compare 
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Figure 9. Comparison with 
for j = {4,5,2} with ( 04 , as, 0,2 

{4>5,4>2) = (16°, 32° 


Ding & Umbanhowar 


(20061 figure 4(b). Phase diagram 
= (3.8,5.8, 0 . 8 )g. [The position of this point for 

is shown in figure [^c).) (a)-(c) stability regions for squares, hexagons 
and quasipatterns respectively. Coloured symbols stable; grey unstable, (d) ‘Most stable’ plan- 


forms: + squares; * hexagons; □ quasipatterns. The r esults are overlaid on f i gure 4 (b) of Ding & 
|Umbanhowar| (2006 1 . Adapted with permission from Ding & Umbanhowar (20061. Copyrighted 
by the American Physical Society. The cross indicates the point which corresponds to the cross 
in figure 


theory and experiment. The experimental results for the pattern at onset are consistent 
with the theory, but sometimes the pattern at onset is not that predicted by considering 
the Lyapunov argument. This is one difficulty in comparing with experiment post festum, 
since in regions of bistability, which pattern is observed can depend on the way that the 
experiments are performed. 

(c) When the additional forcing component (02 ^ 0) is included that promotes the 
mode with frequency 2 d)/ 2 , then in both theory and experiments patterns associated 
with the resonant angle are promoted and consequently appear for a larger region of 
parameter space. 

(d) In the experimental results, the additional forcing causes regions where both hexag¬ 
onal states and super lattice patterns occur close to onset. This is also true in the weakly 
nonlinear analysis, but to a lesser extent, see figure and [m There are two possible 
causes here: firstly, the theoretical results are sensitive to the phase, so changing the 
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Figure 10. Comparison with |Ding & Umbanhowar| (|2006|, figure 2. 
j = {6,7}, (06,0?) = (0,40°),£l)/27r = 16.5 Hz. (a)-(b) stability regions for rectangles 
and hexagons. Coloured symbols s table; grey unstable, (c) ‘Mos t stable’ planform: + rectangles; 
* hexagons, overlaid on figur e 2 of Ding & Umbanhowar (20061. Adapted with permission from 
[Ding fc Umbanhowar| (|2006[) . Copyrighted by the American Physical Society. 
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Figure 11. Comparison with Ding 


j = {6,7,2}, 07 = 7g, (fli6,</>7,02) = (.0,40 


& 


Umbanhowar 
),w/2ty = 


TT 


(|200^, figure 6(a). 
(a)-(c) stability re¬ 


gions for rectangles, hexagons and superlattice patterns respectively. Coloure d symbols stable; 
grey unstable, (d) ‘Most stable’ planforms superimposed on figure 6(a) from Ding & Umban¬ 


howar (2006|) 
American 


Adapted with permission from Ding & Umbanhowar (20061. Copyrighted by the 


lysical Society. The same symbol styles are used for both experiments and theory 
(A superlattice; * hexagons; -I- rectangles; • disordered states). Coloured symbols represent 
theoretical results and black symbols are the experimental results. 


phase can expand the region of stable superlattice patterns. Secondly, the presence of 
the weakly damped 2w/2 mode means that the results are in a region that is close to 
a bicritical point where both 2uj/2 and 4w/2 modes onset simultaneously. In the neigh¬ 
bourhood of this bicritical point, one would expect the regions of stability as predicted 
from a codimension one analysis, at best, be perturbed, and, at worst, be inaccurate. We 
will return to this point in the discussion and note for now the red line in figure that 
marks the onset of this mode. 

(e) Quantitative agreement in most cases only occurs close to onset. 

In addition, in all cases, close to the subharmonic instability boundary the theory 
predicts squares (not shown in the figures), as seen in the experiments. 
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Figure 12. Comparison with Ding & Umbanhowar 
j = {6, 7, 2}, (o6, 07 , 52 ) = (5.2, 7.8,0.6jg, a;/27r = 16.5 


(|2006J), figure 6(b). Phase diagram for 
(a)-(c) stability regions for rectan¬ 


gles, hexagons and superlattice patterns respectively. Col oured symbols stable; grey u nstable, 
(d) ‘Most stable’ planform, superimposed on figure 6(b),|Ding & Umbanhowar (20061. A su¬ 
perlattice; * hexagons; -|- rectangles. Coloured symbols represent theoretical results and black 
symbols are the experimental results. Note that in (d) hexagons are rarely the ‘most stable’ 
planform: whenever hexagons are stable, rectangles usually are too. Although the Lyapunov 
energy for hexagons is often similar to that for rectangles, rectangles nearly always have the 
lowest value. 


6. Discussion 

The Faraday problem is an important pattern-forming system, yet it is particularly 
challenging to tie theory to experiment. The experiments are difficult to perform; the 
parameter regime of interest (large box, moderate viscosity) along with the technical 
difficulties of solving the free boundary Navier-Stokes equations make numerical solution 


of the problem hard, to the extent that although there has been some progress (see Perinet 
et al. (2009, 2012)) it is only very recently that any fully three-dimensional calculations 
of superlattice patterns in the Faraday problem have been reported [Kahouadji et al. 
(2015). The fact that the instabilities result in an entire circle of unstable wavenumbers 


presents considerable theoretical difficulties (see Melbourne (1999)) and has meant that 


theory is, by necessity, restricted to a finite number of modes—the finite number normally 
chosen as the minimal set which will allow for the existence of the experimental pattern 
of most interest. Furthermore, by its nature, weakly nonlinear theory cannot be expected 
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to capture behaviour that is inherently strongly nonlinear: so a weakly nonlinear analysis 
of the transition from non-patterned to patterned state can only be expected to agree 
with experiment close to onset. 


Nevertheless, our detailed comparison with Ding & Umbanhowar (2006) shows very 


good agreement between experiment and weakly nonlinear results at onset, with the 
caveat that the theory often predicts bistability of patterns. Where there is bistability, 
the observed pattern will be dependent on the path taken through parameter space and 
might consist of competition between the different possible stable states. Consequently, 
in some cases it is only possible to say that the weakly nonlinear theory is consistent 
with experiments. 

The results suggest that the qualitative idea that the three-wave resonances determine 
the angles that appear in the patterns does indeed explain many of the patterns observed 
close to onset. 

The results also highlight the particular aspects of the Faraday problem that lead to 
an amplification of bres- Equation (4.5), 


^(^res) — bo 6j., 


brpfi 


G:iG:2 

A? 


suggests that there are two main contributions, one from the quadratic coupling coeffi¬ 
cients («! and 0 : 2 ) and one from the linear damping of the weakly damped mode (A 2 ). 
Promotion of patterns with angle 9res requires the quadratic coefficients to be sufficiently 
large and/or the damping to be sufficiently small. As illustrated in figure [l^a) and (b) 
for j = {4, 5} and in (d) and (e) for j = {4, 5, 2} although increasing 05 and heading 
towards the bicritical point would appear to change the value of the linear damping for 
the 2(i/2 tongue, this is in fact not the dominant effect. This can be seen in figure [T^c) 
and also in (f) where the Floquet multipliers are plotted for both 05 = 0 and a value of 
05 near the bicritical point. It can be seen that increasing the amplitude of the j = 5 
mode in the forcing promotes the 5uj/2 tongue, as expected, but has little impact on the 
damping of the other modes. In contrast, now comparing figure [T^c) with (f), we see 
that the addition of the j = 2 mode does have the effect of reducing the damping of the 
mode associated with the 2w/2 tongue as seen by an increase in the Floquet Multiplier 
from approximately 0.5 for 0,2 = 0 to 0.9 for 02 = 0.8g. 

Rather than changing the damping of the weakly damped mode, increasing the am¬ 
plitude of the j = 5 mode, so heading towards the bicritical point, has the effect of 
increasing the quadratic coefficients ol\ and 0 : 2 : in the limit when 05 = 0 the quadratic 
coefficients are zero. This is because the temporal resonance condition discussed in sec¬ 
tion |4^ is not met when 05 = 0 , so there can be no coupling between the 2 w /2 mode 
and the 4w/2 mode. 


Topaz et al. (2004) consider a weak viscosity limit and show that in the case of 


j = {m, n} forcing the temporal resonance condition results in possible three-wave inter¬ 
actions with an (n — m) mode, where one expects 


brps — 


a\an\ 

IA2I 


( 6 . 1 ) 


where A 2 is the linear damping of the weakly damped mode and a is a coefficient. 
Although this cannot necessarily be expected to apply to the moderate viscosity case of 
the experiments in Ding & Umbanhowar ( |2006 ), we see in figure 14 that the calculated 
values for bres for the j = {4, 5} case have a close to quadratic dependence on 05 as 
suggested by equation (6.1). For j = {m,n,p} frequency forcing that has a three-wave 


resonance with an (n — m) mode, as for the case j = {4,5,2}, Topaz et al. (2004) find 
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Figure 13. Linear stability curves a nd Floquet multipliers for j = {4,5} and j = {4,5,2} 
excitation using parameter values as in |Ding fc Umbanliowar| ( |2006[ ). (a) and (b) Linear stability 
curves for f{i) = a(cosxcos4c5t + sinxcos5t5t), where (a) has x = 0 and (b) has x = 63°. (c) 
the corresponding Floquet multipliers where the solid line is for x — 0 and the dashed line is for 
X = 63°. The Floquet multipliers are calculated for 5 fixed at the value for the minimum of the 
4tj/2 tongue. So, for x = 0, 5 = 4.35g, for x = 63°, 5 = 7A9g (d) and (e) Linear stability curves 
for f{t) = 5(cos xcos 4(1)1 + sin x cos 5ujt + 0,2 cos 2(2)1) where (a) has x = 0 ,02 = aa 2 = O.Sgr and 
(b) has X = 62 °,02 = 002 = 0.8g. (c) The corresponding Floquet multipliers where the solid 
line is for case (d) where 05 = 0 and the dashed line is for (e) where 05 yf 0. The Floquet 
multipliers are calculated for a fixed at the value for the minimum of the 4(2)/2 tongue. So, for 
X = 0, a = 4.40s, for x = 63°, a = 7.27s. 


that the dominant contribution to &res is given by 

^rcs — CrlUnl .^2(^)5 

where ^ = 4>p + 2(j)m — ‘ 24 >n and P 2 is given by 

□ _ IA 2 I + /ra 2 sin$ 


( 6 . 2 ) 


where /r is a coefficient. As 02 —1 0, equation (6.2) reduces to equation (6.1). However, 
for 02 non-zero, the effect of P 2 is to make the leading order dependence of on 02 not 
purely quadratic, so it is not surprising that for this case, an assumption of quadratic 
dependence fits less well. 

Note that the dependence, which comes from a parameter symmetry, does explain 
the strong diagonal structure to the figures showing the pattern dependence as a function 
of two of the phases, see figures and [T^ 

There is a further issue to consider. Increasing the amplitude of the j = 2 mode 
brings the location of the bicritical point where both the 4(Ii/2 and the 2uj/2 mode onset 
simultaneously closer to the pattern onset point. This is seen in figure which shows 
both the linear stability boundaries for both 4(I)/2 and 2ui/2 modes. The proximity of 
the linear stability boundary means that, in order to quantitatively predict the regions 
of quasipatterns, it is likely that one would need to consider the tri-critical problem 
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Figure 14. The value of Kes for j = {4, 5} (crosses) and for j = 14, 5. 21 (dots! for increasing 
as. The solid lines are proportional to a|. Parameter values as in Ding & Umbanhowar (20061 
but with D = 0 . 21 cm^s“^. 


where ACj 12,'bCj 12 and 2Cj/2 modes all occur simultaneously. The onset boundaries for 
the three modes for j = {4, 5, 2} are shown in figure the position of the tri-critical 
point is at approximately ( 04 , 05 , 02 ) = (3.5, 6.4,1.0)^. The situation for j = {6,7,2} 
forcing is similar, and one explanation for discrepancies between the weakly nonlinear 
theory and experimental results may be that they are a consequence of the proximity of 
the codimension three point. This kind of discrepancy was previously demonstrated by 


Riyapan (2012) in a comparison of a codimension one and a codimension two analysis of 


superlattice patterns in a model PDE. 

Another important feature of the analysis considered here is that, the codimension 
one description of the problem results in amplitude equations that are variational. This 
variational structure is useful and it enables us to construct a Lyapunov functional to 
predict the ‘most stable’ pattern. However, this variational framework also means that all 
predicted patterns repeat with either the period or half the period of the drive but cannot 
exhibit more complicated time dependence. As discussed in Rucklidge et al. (2012), in 
the vicinity of the bicritical point and depending on the signs of the quadratic terms, 
an infinite sequence of resonances can take place leading to spatio-temporal chaos. This 
may explain some of the disordered states seen close to onset in the experiments in 


Epstein & Eineberg (2004) and Ding & Umbanhowar (2006). The difficulties of analysing 


and numerically solving the time-dependent Navier-Stokes equations leave a clear role 
for reduced models here. For example, the model of Zhang & Vihals (1997) is derived 
from the Navier-Stokes equations under the assumptions of infinite depth and weak 
viscosity but, as was shown in Skeldon & Porter (2011), weakly nonlinear analysis of 
the Zhang-Vinals equations agrees surprisingly well with weakly nonlinear analysis of 



























26 


A. C. Skeldon and A. M. Rucklidge 



aiig) 0 0 


Figure 15. Surfaces for the onset of modes with frequ encies AujI2,bu} 12 and 2dj/2. Parameter 


values as in 


Ding & Umbanhowar (20061 but with v — 0.21cm^s 


the full Navier-Stokes equations even for the moderate viscosity values used in many 
Faraday experiments. 

Finally, our results suggest that weakly nonlinear theory is a useful tool to understand 
the pattern transitions near onset, but there are some caveats. Ideally experiments and 
theory would be done hand-in-hand otherwise it may be hard to get agreement even in 
the linear theory, as we saw in section 3; or establish which modes are present in patterns, 
as we saw in section 4; or establish if there is bistability, as we saw in section 5. There are 
even problems at the basic level of the definition of the forcing in equation (2.1). In single 
frequency experiments changing the sign of f{t) is equivalent to a phase shift in time, so 
the sign of the forcing is not important. However, for multiple frequencies the analogous 
result is perhaps less transparent: changing the sign requires not just a translation in 
time but requires an altered definition for the phases. The particular sign convention 
used by any one experimental group will depend on how the accelerometer is wired up to 
the experiment. Specifically, it will depend on whether or not the accelerometer records 
its maximum value when the oscillator is at its highest point or at its lowest. This level 
of detail is not normally recorded — perhaps because an assumption is made that it does 
not matter. However, this can lead to some confusion: for example, the phase invariant $ 
discussed in Porter et al. (2004| is derived with one particular choice of sign for f{t) and 
thus one specific definition of the phases for the forcing components. [Ding fc Umbanhowar] 
(2006) note that they find a difference between the position of the maximum value of the 
phase invariant and the theory of Porter et al. (2004) and between their results and those 
of Epstein & Fineberg (2006), but in both cases these differences can be explained by a 


different convention for the sign of the forcing f{t). Based on the observations in Ding 


& Umbanhowar (2006), we have used the same sign convention as Porter et al. (2004) 
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for our comparison with the work of 

Epstein & Fineberg 

(2006 

) and the opposite sign 

convention for Ding & Umbanhowar 

2006 

)■ 



Appendix A. Calculating Floquet multipliers 

The Floquet multipliers can be calculated by considering the linear stability of the triv¬ 
ial solution (2.9) of equation (2.6) and associated boundary conditions to perturbations 
of the divergence free form 


This leads to a fourth order equation for W{z,t) for —h/l < z < 0 

{dt - C {-e + 5 ,,)) {-e + d,,) w = o, 

with boundary conditions at z = —h/l, 

Wi-h/l,t) = d,WU^j = 0 


(Al) 


and at z = 0, 


dtZ = W, 

(fc2 + 5,,) W = 0, 

(dt + C {3^ - d,,)) d,W = - ((1 + f{t))k^ + Bk^) Z. 


A finite difference discretization of equation ( |A 1[ ) and its boundary conditions is then 
carried out by letting W" = W{zj,tn),j = O..J,Zj = —h/l + hj/lJ and = n6t and 
Z"^ = Z{tn). The resulting map is of the general form 



/ \ 


/ wr \ 




VF 2 " 

A 


= B(t„) 



W/+1 


wy_i 


M/J+i 


wy 


\ / 


\ / 


(A 2) 


where A and B are (J -I- 1) x 


(J -I- 1) matrices given by 


/ r 1 0 

1 r 1 

0 1 r 


A = 


0 


0 

1 

1 

0 


\ 


0 


r 1 0 

2 r 0 

0 1 
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and 


f c + C5t/5z^ 
b 

Cbtjbz^ 

0 


B = 


V 


where 


b 

c 

b 

CStlSz'^ 


CSt/Sz 


0 

CStjSz^ 

b 

c 


0 C6t/5z‘^ 


CSt/5z 


CSt/5z 


0 

CStlSz^ 


c 

b 

2CStlSz^ 


c — Cbtjbz^ 
6 + e 
0 


r = — (2 + k^bz^) 
6=1- 2Ck'^bt - AC 


hk^bz^) + 
(2 - k^bz‘ 


b^ 
bz^ 
bt 


i& + Ak^bz^ + k^bz'^) 
dz^ ^ ^ 


c = - (2- 

e = 1 — ACk^bt 

r = "^btbz ((1 + /(i„)) + Bk^) 

9 = -cA(2-fc»fc»)^ 

The matrix B contains time dependence through the term /"(tn)- 

The map (A 2) is iterated through one period, T, of the drive resulting in a map 
= DW° that takes W° to where N = T/bt. The eigenvalues of the matrix D 
then give the required Floquet multipliers. 


This paper builds on discussions about pattern formation in the Faraday problem with 
many individuals over the years. In particular the authors would like to acknowledge the 
input of Mary Silber, Jeff Porter and Jay Fineberg. We thank the referees for their 
thorough reading of the paper and their helpful comments. 
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